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On" 

' The near membrane layer is defined as a region where the concentration of the substance trans- 

^\ , ported across the membrane drops k times. The time evolution of such a layer is studied experi- 

mentally by means of the laser interferometric method. It is shown that within the experimental 
errors the thickness of the near membrane layer grows in time as with the coefficient a being 
, independent of the initial concentration and the membrane permeability. Time evolution of the near 

' membrane layers is also analyzed within the theoretical approach earher developed by one of us. 

The regularities found experimentally fully agree with the theoretical expectations. 
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O ■ I. INTRODUCTION 

oo 

' The transport in membrane systems is of great interest in several fields of technology ||l[ , where the membranes 

■ are used as filters, and biophysics ^, where the membrane transport plays a crucial role in the cell physiology. 
The diffusion in a membrane system is also interesting by itself as a nontrivial stochastic problem. While the time 
dependent concentration profiles of the substance transported across the membrane give a detailed description of the 

'"O ' macroscopic substance motion, we are often interested only in the regions with a sufficiently large concentration. Such 
, a situation occurs when the phenomenon under consideration strongly depends on the concentration. For example, 
^ ' the hydrodynamic stability in the membrane systems studied in |^ is controlled by the Rayleigh number, which in 
turn depends on the transported substance concentration |Q. Therefore, it is sometimes convenient to introduce the 
so-called near membrane layer (NML) where the substance concentration drops k times ||^. In this paper we study 
NML experimentally and theoretically. 

Our experimental investigation is carried out by means the laser interferometric method. The laser light is spilt into 

■ two beams. The first one goes through the membrane system parallelly to the membrane surface while the second, 
reference one goes directly to the light detecting system. The interferograms, which appear due to the interference 
of the two beams, are controlled by the refraction coefficient of the solute which in turn depends on the substance 
concentration. The analysis of the interferograms allows one to reconstruct the time dependent concentration profiles 
of the substance transported across the membrane. Further one can find how NML evolve in time. Wc show that the 
time evolution of the NML thickness manifests surprisingly simple regularities. 

The time dependence can be studied using the Smoluchowski (diffusion) equation. However, one has to impose two 
boundary conditions at the membrane surface. The first one is provided by the substance current conservation but 
there is no obvious choice of the second condition. When the membrane has a finite thickness, the diffusive transport 
within the membrane is often described by the Smoluchowski equation as well Q. The diffusion constant however 
differs from that one which is in the regions outside the membrane. Then, one assumes B that the ratio of the 
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concentrations at both sides of each of the membrane surfaces equals a constant which is a free parameter. We find 
this approach as not very satisfactory. Since a real membrane is not homogeneous and its internal structure is rather 
complicated, using the diffusion equation inside the membrane is rather questionable. The boundary condition, which 
fixes the concentration ratios, is introduced without physical justification. When the membrane thickness goes to zero, 
the membrane selectivity vanishes entirely, and consequently, the approach is hardly useful for very thin membranes. 
In the series of papers one of us has developed the approach which is also based on the Smoluchowski equation. 

The boundary condition is well motivated and the membrane can be treated as an infinitely thin wall characterized 
by the permeability coefhcient. The time dependence of NML can be easily derived and, as we show here, it agrees 
very well with the experimental data. 

The paper is organized as follows. In Sec. II we present the experimental procedure with the results on the 
concentration profiles and the near membrane layers. Sec. Ill is devoted to the theoretical considerations. The 
solution of the Smoluchowski equation is found and the time dependence of the NML thickness is derived. The 
predictions of our theoretical model are compared with the experimental data. We summarize our study in Sec. IV. 



The membrane system under study is the cuvette of two chambers separated by the horizontally located membrane. 
Initially we fill the upper (lower) chamber with the aqueous solution of the ethanol while in the lower (upper) one 
there is pure water. Then, the ethanol diffuses to the lower (upper) chamber. Since the concentration gradients are 
in the vertical direction only, the diffusion is expected to be one dimensional (along the axis x). In other words, the 
ethanol concentrations are assumed to be uniform in the planes parallel to the membrane. 

As already mentioned we employ the laser interferometric method to measure the time dependent concentration 
profiles in the membrane system. Let us note that the measurement does not disturb the system under study. The 
experimental setup is described in 0, here we only mention that it consists of the measurement cuvette with the 
membrane, the Mach-Zehnder interferometer |l^ ] including the He-Ne laser, TV-CCD camera, and the computerized 
data acquisition system. 

The interferograms are sensitive to the variation of the refraction coefficient within the membrane system and 
consequently to the concentration gradients. When the solute is uniform the interference fringes are straight and they 
bent when the concentration gradient appears. The example of the interference images are shown in Fig. 1. The 
substance concentration at x is determined by the deviation d of the fringes from their straight line run. Since the 
relation between the concentration C and the refraction coefficient is assumed to be linear, we have 



where Cq is the initial substance concentration; a is the proportionality constant between the concentration and the 
refraction index, a = 3.19 • 10''^ mo\/m^ for the ethanol aqueous solution; A is wavelength of the laser light; h denotes 
the distance between the fringes in the field where they are straight lines; / is the thickness of the solution layer in the 
measurement cuvette. Recording the interferograms with a given time step one can reconstruct the time dependent 
concentration profiles. 

We measured the profiles for several values of the initial ethanol concentration and for two cellulose membranes of 
different permeabilities. The membrane thickness is, respectively, 0.17 mm and 0.01 mm. In Fig. 2 we present the 
concentration profiles taken at different moments of time. Repeating several times the measurements at the same 
conditions we tasted the stability of our results. While the shape of the concentration profiles has appeared to be 
rather stable, the absolute normalization has varied by about 15 %. The error bars shown in Fig. 2 just correspond 
to this uncertainty. 

Since the specific weight of the ethanol is significantly lower than that of water, there can appear a convective 
motion in the system when the ethanol is initially in the lower chamber. To check whether the gravitational force 
significantly influences the process we have studied the ethanol transport from the upper chamber to the lower one 
and then the opposite configuration. We have found that for sufficiently low initial ethanol concentrations under 
study (125 - 750 mol/m'^), the two configurations are within the experimental uncertainties equivalent to each other. 

Having the profiles one can define the near membrane layer (NML). When the substance diffuses across the mem- 
brane into the pure solvent, the thickness I of NML is defined as a length at which the concentration decreases k times 
i.e. 



II. EXPERIMENT 
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with a; = being the membrane position]^ 

Taking k = 12.5 we found the thickness of NML from the earher obtained concentration profiles. In Figs. 3 and 4 
we present the thickness as a function of time. As seen it manifests the remarkable properties. The NML thickness 
appear to be independent, within the experimental errors, of the initial concentration (Fig. 3) and of the membrane 
permeability (Fig. 4). The thickness grows in time as As we have already mentioned our concentration profile 
measurements suffer from the absolute normalization uncertainty of about 15 %. However, this uncertainty does not 
influence the NML thickness because the absolute normalization coefficient drops out entirely in the NML defintion 
(|l|). Therefore, the estimated errors of the data points shown in Figs. 3 and 4 are contained within the data point 
symbols. In the next section we discuss the NML thickness from the theoretical point of view. 



III. THEORY 



Let us briefly present the approach to the membrane transport developed in [^-10|. The concentration profile 
C{t,x) is assumed to satisfy the Smoluchowski equation 



dt dx^ ' 



(2) 



where D is the diffusion constant. As already mentioned we have experimentally checked that the we deal with a pure 
diffusion. Therefore, in contrast to our earlier study ||], the term responsible for the convection is neglected in (||). 
The membrane is treated as an infinitely thin, partially permeable wall. Then, one needs two boundary conditions at 
the wall to solve eq. (|^) . The first one is provided by the substance current conservation. The second one is given by 
the relation 



J{t,x = Q) = {l-5) J°(t,a: = 0) , 



(3) 



where Jit, x) — —D dC jdx is the substance current in the membrane system while J°(<, a;) = —D dC^ /dx denotes the 
current in the identical system but with removed membrane; 5 is the membrane permeability coefficient, < (5 < 1. 

As already mentioned, the solute concentration is initially zero in one chamber of the membrane system and it is 
finite and uniform in the other one i.e. 



C{t = 0,a;) 



Co 




for a: < 
for a: > 



The solution of eq. (0) for the given boundary and initial conditions reads 



C{t,x)^Co 



1 



■ erfc ( — 



for a; < 



and 



C(i,a;) =Co^erfc(- 



for a; > 



2 y2^/Dt' 

with erfc(a;) being the complementary error function defined as 

2 f°° _ 2 

erfc(a;) ~ —= / dte ^ . 
Vi" Jx 

Substituting the solution (||) into the layer thickness definition (|^) we get 

l{t) =aVi , 



(4) 



(5) 



(6) 



'^If the substance diffuses to the region of the nonzero initial concentration Co, the definition (0) is generalized as C{t, 0) — Co 
k[C{t,l) - Co]. This form can be applied to the layer on both sides of the membrane. 
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where the coefRcient a depends solely on D and k as 



a = 2/Derfc"^fi 



and erfc"^(0.08) ^ 1.24. One sees that the relation (^) fully agrees with the experimental results shown in Figs. 3 
and 4. Indeed, the dependence on the initial concentration Co and the membrane permeability coefficient S, which is 
present in the solution (^), drops out entirely in eq. (||). 

The near membrane thickness is independent of Co due to the linearity of the diffusion equation(^). The 
cancellation of the membrane permeability in eq. (^) is much less trivial. A solution of the diffusion equation (|^) 
usually depends on / Dt, which is the only dimensionless combination of D, x and t. Then, any length must be 
proportional to ^/Dt as in our eq. (|^). However, the diffusion coefficient D is not the only one dimensional parameter 
in the problem under consideration. For example, the membrane thickness d can play such a role. Then, we can 
construct several dimensionless combinations of x and t e.g. xd/ Dt. In our model the membrane is infinitely thin, so d 
is irrelevant. However, the dimensional parameter can be introduced to the problem through the boundary condition. 
For example, the dimensional membrane permeability parameter x appears in the condition of the form 

J{t,x = 0) = x\c{t,x^ 0+) - C(i,x = 0") 

Then, the solution of the diffusion equation, which satisfies this boundary condition, depends not only on x^ jDt but 
on other dimensionless combinations such as xV^- Then, the thickness of the near membrane layer is x dependent. 
Thus, one sees that our formula (^) is independent of the membrane permeability because of the specific choice of the 
boundary condition (|^). 

In our previous studies we took the diffusion coefficient form the literature to fit the data. However, the 
coefficient is known to be sensitive to the temperature and the ethanol concentration in water. One sees that our 
formula (||) solely depends on the diffusion coefficient. Therefore, D can be obtained directly from our data (shown 
in Figs. 3 and 4) by means of a single parameter fit. In this way we get D = 0.76 • 10^^ m^/s. Since the value 
of D is known, there is only one parameter 5 to be fitted when the solution (||) is compared with the experimental 
concentration profiles. The resulting theoretical curves are shown in Fig. 2. The agreement with the experimental 
data is reasonably good but a comment is in order here. 

Our theoretical approach assumes that the membrane is infinitely thin while the membranes which have been used 
in the measurements are, obviously, of the finite thickness. Further, there is a near membrane dead zone where the 
concentration measurement is unreliable or even impossible. The dead zone appears due to the imperfection of the 
cuvette edge and the small deformations of the membrane during the measurements. The latter effect, which is much 
more important for the thinner membrane, is not very well controlled in our set-up and it leads the substance stirring 
in the very vicinity of the membrane surface. To take into account all these effects we have shifted the theoretical 
curves from Fig. 2 by 0.3 mm to the left. This is, of course, not a quite satisfactory procedure and in the future 
studies we intend to resolve the problem in a systematic way. The main objective of this study, however, are not the 
concentration profiles but the near membrane layers. We have carefully checked that the time dependence of the near 
membrane layer is infiuenced insignificantly when the solution (|^) is shifted by 0.3 mm. 

IV. FINAL REMARKS AND SUMMARY 

The solution of the diffusion equation found with the boundary condition (^) is identical with the solution which 
satisfies the boundary condition 



C(i,x = 0") = KC(t,a; = 0+) (7) 



when K is related to 8 as 



Therefore, the formula can be obtained within the boundary condition which is well known in the diffusion 
theory. However, as far as we know the time dependence has not been earlier discussed. 

We also mention that the boundary condition applies for a pure diffusion i.e. when the convection term is 
neglected in eq. (||). However, when the convection is present in the system, one easily generalizes eq. (||) [^-p^] 
and introduce the parameters which separately determine, as in the Kedem-Katchalsky approach the membrane 
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permeability with respect to the diffusion and to the convection. On the other hand, we see no way to generalize eq. 
(0) to distinguish between the diffusive and convective membrane permeability. Therefore, we would like to stress 
that the equivalence of eqs. (|^) and holds only for a pure diffusion regime. 

Let us summarize our study. Using the interferometric method the concentration profiles have been measured 
for the aqueous solution of the ethanol diffusing across the membrane into pure water. The time evolution of the 
near membrane layer has been then analyzed. The thickness of the layer, which is defined as a length where the 
concentration drops k times, appears to grow in time as a\/i with the proportionality coefficient a being remarkably 
independent of the initial ethanol concentration and of the membrane permeability. While the independence of Cq is 
the result of the problem linearity, the origin of the cancellation of the membrane permeability is much less trivial. 
Theoretical analysis performed within the approach proposed by one of us has shown that a depends only on k and 
the diffusion constant D due to the specific choice of the boundary condition at the membrane. We have used this 
fact to obtain the numerical value of D fitting the experimental time dependence of the near membrane layer. Then, 
the experimentally found concentration profiles have been fitted by the solution of the Smoluchowski equation. The 
membrane permeability coefficient S is then a single fit parameter. 
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Figure Captions 



Fig. 1. The interferograms which are analyzed to obtain the concentration profiles. There is initially uniform ethanol 
solution of the concentration 125 mol/m^ in the upper part of the measurement cuvette. The interferograms are taken 
at several values of time: (a) - 240 s, (b) - 600 s, and (c) - 1200 s. 

Fig. 2. The concentration profiles for Co = 125 mol/m^ taken at several values of time: o - 240 s, A - 600 s, and 
□ - 1200 s. The permeability of the membrane from (a) is significantly smaller than that from (b). The solid lines 
represent the Smoluchowski equation solution with D = 0.76 • 10~^ m^/s and 6 = 0.51 for (a) and 6 = 0.11 for (b). 
The theoretical curves are shifted to the left in (a) and (b) by 0.3 mm (see text). 

Fig. 3. Time evolution of the near membrane layer for several values of the initial ethanol concentrations: o - 
125 mol/m^, A - 250 mol/m^, □ - 500 mol/m^, and ★ - 750 mol/m^. The solid line represents l{t) = a\/i with 
a = 6.8 -10-5 m/sV2. 

Fig. 4. Time evolution of the near membrane layer for two membranes of different permeabilities: o - membrane 1, 
A - membrane 2. The solid line represents l{t) = a\ft with a = 6.8 • 10~^ m/s^/^. 
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